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The  theoretical  background  for  the  study  of  hydro- 
elastic  problems,  recently  Introduced  to  the  designers  of 
naval  craft,  has  been  taken  mostly  from  the  analogous  field 
of  aeroelastlclty.  Two  problems  In  hydrofoil  design  that 
have  received  some  attention  are  flutter  (oscillatory  diverg¬ 
ent  motion)  and  divergence  (exponentially  divergent  motion). 

Tlie  boundary  be  I. ween  sizable  and  unstable  motion  In  each  case 
Is  the  critical  flutter-speed  or  cid  tlcal  divergence-speed. 

This  volume  discusses  only  critical  flutter-speed. 
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ABSTRACT 


The  theoretical  procedures  commonly  used  by  aero- 
elastlclans  were  applied  to  predict  the  flutter  speed  of  a 
rigid  hydrofoil  that  had  two  degrees  of  freedom.  The  results, 
compared  with  corresponding  experimental  measurements,  Indi¬ 
cated  a  discrepancy  between  theoretical  and  experimental 
flutter  speeds  at  low  density  ratios;  the  predicted  asymptotic 
behavior  of  flutter  speeds  occurred,  but  at  a  lower  density- 
ratio.  In  addition,  the  accuracy  of  the  circulation  terms  Is 
more  doubtful  than  that  of  the  added  mass  and  linear  terms  In 
the  theory. 
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I. 


INTRODUCTION 


Destructively  large  stresses  or  undesirable  levels 
of  vibration  have  been  induced  In  airfoils  operating  near  the 
critical  flutter-speed.  Theoretical  and  experimental  Investi¬ 
gations  have  shown  that  for  certain  elementary  elastic  con¬ 
figurations,  flutter  of  foils  In  water  Is  unlikely  (references 
1,  2,  and  .  Figure  1  shows  the  range  of  density  ratios 
applicable  to  these  simple  hydrofoil-configurations.  The 
magnitude  of  the  dynamic  pressure  encountered  by  submerged 
control  3urfaces--rudders,  stabilizing  -  fins,  hydrofoils,  bow 
planes,  etc, --is  greater  than  that  encountered  by  airfoils, 
except  for  those  used  on  recent,  high-speed  aircraft.  There¬ 
fore,  one  Is  led  to  the  following  question:  why  have  hydro- 
dynamic  control  surfaces  not  experienced  flutter? 

In  particular,  bending -torsion  flutter  of  canti¬ 
lever-supported  hydrofoils  has  been  investigated  (reference  3) • 
In  this  Investigation,  a  critical  density-ratio  was  predicted 
below  which  flutter  was  not  possible.  Isolated  experiments 
using  simple  supports,  confirmed  these  results.  In  these 
experiments,  no  flutter  was  obtained  and  none  predicted. 

However,  the  low,  overall  damping  associated  with  the 
unsteady  hydrodynamilc  forces  on  the  rudders  of  a  destroyer, 
which  recently  experienced  severe  and  sustained  hull  vibrations, 
has  provided  the  Impetus  for  further  Investigations  of  hydro- 
dynamic  flutter  (reference  5)-  Cursory  Investigations  of  this 
vibration  Indicate  that  more  realistic  support  conditions,  and 
other  conditions,  may  reduce  the  lowest  value  of  density  ratio 
at  which  hydrofoils  will  flutter.  Therefore,  the  reliability 
of  theoretical  predictions  of  flutter  speed  at  low  values  of 
density  ratio  must  be  determined. 
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In  previous  hydroelastic  work.  It  was  assumed  that 
aerodynamic  theories  would  apply  directly  to  the  hydrodynamic 
problem.  Assuming  this  to  be  true,  these  theories  should  cor¬ 
rectly  predict  the  flutter  speed  of  a  hydrofoil.  This  Investi¬ 
gation  was  conducted  to  experimentally  verify  the  accuracy  of 
aeronautical  theories  when  used  to  predict  flutter  speed 
(throughout  the  range  of  parametric  values  of  Interest  In 
hydrofoil  design)  . 

Some  previous  work  along  these  lines  has  been  per¬ 
formed.  The  NACA  Investigated  flutter  of  light, cantilever- 
supported  wings  In  an  alrstream  made  heavy  by  the  use  of  Freon 
and  air  mixtures  (reference  ^l)  .  Tlie  experimental  points  In 
Figure  1  are  representative  of  the  results  obtained  by  NACA. 

The  range  of  density  ratio  In  this  Investigation  did  not  ex¬ 
tend  into  that  for  a  cantilever-supported  hydrofoil.  However, 
the  range  ol'  density  i-atlo  In  wlilch  flutter  was  obtained  ex¬ 
perimentally  extends,  l^y  a  small  amount  below  the  critical 
density-ratio,  into  the  range  in  which  flutter  Is  theoretical¬ 
ly  impossible.  It  Is  this  discrepancy  which  was  Investigated. 

In  view  of. the  inconclusiveness  of  previous  experi¬ 
mental  and  theoretical  results,  a  single  set  of  experiments 
was  conducted  over  a  range  of  density  ratio  extending  from 
region  1  through  regions  2  and  3  and  well  into  region  4 
( Figure  l) . 

In  these  expei'iments,  speed  and  density  ratio  were 
varied;  all  other  parameters  were  constant.  Another  series 
of  experiments  Is  being  conducted  at  several  additional 
values  of  center- of-gravlty  location  and  radius  of  gyration 
to  Investigate  further  the  discrepancy  between  theory  and 
experiment . 


II. 


TEST  SETUP 


A.  MODEL  AND  SUPPORTS 

The  experimental  setup  was  designed  so  that  flutter 
would  be  obtained  under  controlled  conditions  that  closely 
duplicated  the  theoretical  assumptions;  the  setup  did  not  typi¬ 
fy  any  hydrofoil  application  (Figures  2  through  6) .  The  ex¬ 
periments  were  performed  In  the  High-Speed  Facility. 

The  model  had  a  chord  length  ( c)  of  6  Inches  and  a 
span  of  12  Inches.  The  profile  was  a  thin  symmetrical  shape, 
(NACA  0012) .  The  offsets  for  the  model  were  obtained  from 
reference  6.  The  material  used  for  the  model  was  a  plastic 
whose  properties  are  summarized  In  Table  I. 

1 .  End  Plates 

End  plates  were  used  to  provide  two-dimensional  flow 
and  were  attached  to  the  carriage  by  struts  and  aluminum  box- 
beams.  The  chord  plane  of  the  model  was  arranged  vertically 
and  the  end  plates  were  above  and  below  the  model,  each  with  a 
maximum  clearance  oC  0.006c.  The  end  plates  were  11.29c  long 
by  6.71c  wide.  The  rotational  axis  of  the  model  was  located 
In  the  middle  and  ^l.31c  aft  of  the  leading  edge  of  the  end 
plates.  The  end  plates  were  made  of  l/8-lnch  aluminum  plates 
stiffened  In  both  the  transverse  and  longitudinal  direction. 

The  transverse  stiffening  members  were  covered  with  aluminum 
sheets,  which  provided  fairing  to  reduce  drag. 

2 .  Sting  Support 

The  model  was  supported  by  a  vertical  sting  that 
passed  through  the  upper  end-plate  and  the  water  surface  to  a 
flexure  balance  that  was  entirely  above  water.  The  sting  was 
a  stainless-steel  tube  (l-lnch  OD,  0.87-inch  ID,  DWG  No.  I6) 
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protected  from  the  stream  velocity  In  the  region  between  the 
upper  end-plate  and  the  water  surface  by  a  faired  shield. 

To  reduce  the  effect  of  the  hole  In  the  upper  end- 
platej  a  small,  circular  plate  made  of  l/l6-lnch  stainless 
steel  was  attached  to  the  lower  end  of  the  sting.  The  radius 
of  this  auxiliary  end-plate  overlapped  the  hole  In  the  upper 
end-plate  when  the  sting  was  displaced  to  either  side.  The 
auxiliary  end-plate  was  below  the  main  end-plate  and  had  a 
maximum  clearance  of  0.006c. 

3 .  Balance 

The  sting  was  attached  to  the  fixed  support  through 
a  flexure  balance  (Figures  5,  6,  and  7)  •  Two  modes  of  motion 
for  the  model  in  the  horizontal  plane  were  permitted  by  this 
balance;  rotation  about  a  vertical ,  spanwlse  axis  and  horizontal 
translation  nonnal  to  the  direction  of  motion.  Mechanical 
stops  limited  these  motions  to  ±3/^l-lnch  of  translation  and 
-2  degrees  of  rotation.  The  axis  of  the  sting  was  In  line 
with  the  quarter-chord  axis  of  the  model  and  was  located  at 
the  apex  of  a  pair  of  V- frames  (Figure  7) .  These  frames, 
located  In  parallel  horizontal  planes  one  above  the  other, 
provided  a  linear  elastic  restoring  moment  In  the  desired 
rotational  degree  of  freedom  and  were  effectively  rigid  In  all 
other  degrees  of  freedom.  These  frames  were  I-beams  with 
lightening  holes  In  the  webs  and  were  necked  down  at  both  ends 
as  shown  In  Figure  7.  Most  of  the  flexibility  of  these  beams 
was  concentrated  In  the  necked-down  sections  (flexures);  the 
dimensions  of  the  flexures  were  chosen  to  give  the  desired 
rotational  stiffness.  The  base  of  the  V-frames  was  connected 
to  the  fixed  support  through  four  parallel  bars,  which  allowed 
the  whole  unit  to  translate  In  the  desired  direction.  In  a 
manner  similar  to  the  members  of  the  V-frames,  the  ends  of  the 
parallel  bars  were  necked  down  to  provide  the  desired  trans¬ 
lation  flexibility  (Figure  '()  .  The  translation  bars  were 
aluminum  channels  and  the  rotation  beams  were  machined  alumi¬ 
num  stock. 
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A  spring-loaded,  locking  mechanism  was  Installed 
that  held  the  model  In  the  desired  Initial  position  until  It 
was  released  by  a  solenoid  controlled  by  a  switch  at  the 
operators  station. 

The  flexure  balance  was  anchored  to  the  carriage  by 
a  supporting  structure.  To  provide  greater  rigidity  against 
rotation  of  the  apparatus  In  the  vertical  transverse  plane, 
an  outrigger  was  attached  to  an  auxiliary  wheel  that  ran 
along  the  side  of  the  tank. 

The  apparatus _ satisfied  all  requirements  for  the 
support  of  the  model.  The  rotation  spring  stiffness  came  with¬ 
in  8^  of  the  desired  value  and  the  translation  spring  stiffness 
came  within  3^’-  These  errors  arose  mainly  from  the  estimate 
of  the  effect  of  fillets  at  both  ends  of  the  flexures.  A  flex¬ 
ure  length  equal  to  the  actual  I'lexure  length.  Including  both 
fillets,  minus  one  fillet  radius,  gave  values  for  predicted 
stiffnesses  within  2%  of  the  measured  values.  The  rigidity  of 
support  against  undeslred  motions  was  found  to  be  well  within 
reasonable  limits.  The  frequency  of  vertical  translatory 
oscillations  of  the  model  in  air  was  found  to  be  three  times 
greater  than  the  highest  frequency  encountered  during  the 
tests  In  water.  All  other  observed  I'requencles  of  extraneous 
motions  were  five  or  more  times  the  highest  test  frequency. 

B.  INSTRUMENTATION 

The  motions  of  the  model  were  measured  by  recording 
the  unbalance  of  strain-gage  bridges.  Strain  gages  were  mount¬ 
ed  on  the  flexure  balance  and  the  signal  outputs  from  the 
strain-gage  bridges  were  recorded.  To  measure  the  speed  of 
the  apparatus,  an  electronic  counter  was  used  to  determine  the 
time  required  for  the  apparatus  to  pass  between  two  locations 
a  known  distance  apart.  In  addition,  a  tachometer  generator, 
attached  to  the  carriage,  was  activated  by  a  wheel  running  on 
the  carriage  rail  to  measure  the  Instantaneous  speed  of  the 
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apparatus.  The  average  speeds  obtained  by  these  two  methods 
were  in  agreement  and  no  appreciable  speed  fluctuations  dur¬ 
ing  the  runs  were  observed. 

In  general,  slide  rule  accuracy  was  maintained 
throughout  the  experimental  analysis.  Scatter  in  the  data 
is  a  result  of  actual  extraneous  effects  in  the  experiments 
and  Is  not  caused  by  Inaccuracies  In  the  Instrumentation. 

The  parametric  values  obtained  In  the  experiments 
are  shown  in  Table  II.  In  view  of  the  large  amount  of  weight 
added  externally  and  rigidly  to  the  foil,  the  results  of  the 
experiments  described  here  should  not  be  viewed  as  represen¬ 
tative  of  any  hydrofoil  application.  In  practice,  density 
ratios  greater  than  1.0  are  seldom  encountered,  If  ever. 

The  experimental  flutter  speeds  and  frequencies  are 
given  In  Table  II  and  the  results  are  plotted  In  Figures  10 
■through  2>l.  The  corresponding  theoretical  results  are  pre¬ 
sented  along  with  these  results. 

C.  l^IGHTS 

Two  sets  of  weights  were  connected  rigidly  to  the 
sting  to  obtain  the  selected  center-of-gravlty  (c.g.)  loca¬ 
tion  radlus-of-gyratlon  (r^b),  uncoupled  natural  fre¬ 

quency  ratio  (oD  /o)  ),  and  to  vary  the  density  ratlo(jLi  )  while 

1 1  w  1 1 

maintaining  all  other  parameters  as  constants  (Table  II) . 

The  uncoupled  natural  frequency  In  each  mode  of 
vibration  can  be  defined  as  follows: 

=  ti^/m  and  coj  =  \  l) 

where  =  support  stiffness  In  translation  per  unit  span 
=  support  stiffness  In  rotation  per  unit  span 
m  =  total  oscillating  mass  per  unit  span 

I^  =  mass  polar  moment  of  inertia  of  rotating  parts 

about  the  rotational  axis,  per  unit  span. 
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The  dlmenelonless  radius  of  gyration  can  be  defined 
as  follows: 

(2) 

where  b  Is  the  half-chord  length  of  the  model.  When  equa¬ 
tions  1  and  2  are  combined,  the  uncoupled  natural  frequency 
ratio  Is 

(“lA)"  =  (V\'  (''a'’)"  (3) 

Thus,  the  stiffness  ratio  (Kj^/K^)  Is  determined  by  the  values 

of  cc./jo  ,  r^,  and  b. 
h  a  a 

Let  the  subscript  b  denote  the  characteristics  of 
the  apparatus  without  weights  and  the  subscript  a  denote  the 
characteristics  of  the  apparatus  with  the  first  set  of  weights. 
Let  the  subscript  n  =  0,  1,  2,  ...  denote  the  desired  values 
of  the  properties  numbered  successively,  starting  with  the 
smallest  mass  condition.  Therefore,  the  necessary  additional 
mass  (m^)  can  be  determined  as  follows: 

■"a  “  %  ■  ”b  = 

where  =  density  ratio  of  Initial  mass-condition  =  m^Trpb^ 
p  =  mass  density  of  water. 

To  find  the  required  location  of  the  c.g.  of  m  ,the  moment 
of  mass  about  the  selected  c.g.  location  Is  set  equal  to  zero 
as  follows: 

^a™a  ^b'"b  '  °  (5) 

where  x  Is  the  distance  from  the  selected  c.g.  location  to 
the  c.g.  of  the  mass  Indicated. 

Equations  4  and  5  were  used  to  determine  the  amount 
of  mass  necessary  to  obtain  the  Initial  mass-condition  and 
the  location  of  the  c.g.  of  this  mass.  The  distribution  of 
m  was  selected  to  provide  the  required  mass  polar  moment  of 
Inertia  Therefore,  two  weights,  each  with  mass  m /2, 

were  symmetrically  placed  on  either  side  of  the  c.g.  location 
(x  )  .  The  distance  between  the  masses  was  adjusted  until 
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OC  Q  0l>  o 

In  this  manner,  the  desired  values  of  x^,  r^, 

In  addition  to  the  Initial  value  of  density  ratio  ( Mq) ^  were 
obtained.  The  second  set  of  weights  were  designed  to  provide 
the  selected  values  of  ju^  while  maintaining  a,  x^,  r^,  and 
constant. 

The  position  of  the  rotational  axis  ( ab)  and  the 
stiffness  ratio  were  fixed;  these  properties  remained 

unchanged  throughout  the  experiments.  The  desired  x^  was 
obtained  with  the  first  set  of  weights.  To  keep  the  c.g. 
location  constant,  the  c.g.  of  the  second  set  of  weights  was 
put  at  the  selected  c.g.  location  and  thus  did  not  change 
x^.  Equation  3  shows  that  will  remain  constant  as  m 

Increases  when  r  Is  constant.  Therefore,  the  ratio  I„/m 

'Jj  vv 

was  kept  constant  as  the  weights  were  added. 


At  each  mass  condition  (n),  the  moment  of  Inertia 
(loj^)  and  the  mass  (m  )  may  be  divided  Into  the  sum  of  the 
respective  values  i'or  the  lowest  mass  condition  plus  the 
Increment  due  to  each  weight  added  (l  =  1,  2,  ...,  n  -  l). 


n 


'a 


-  I,  +  E 

a  ^ 


n 


1-1 


a, 


n 

m  =  m  +  E 
^  °  1=1 


Am, 


n  =  1,  2.  .  . 


(7) 


By  substituting  these  quantities  Into  the  definition  of  the 
radius  of  gyration,  equation  2  yields 


n  n 

I  +  E  AI  =  (r  b)^  Cm  -I-  E  Am.  1, 
“o  1=1  ^1  °  1=1  ^ 

and  subtracting  equation  6  from  equation  8  leaves 

Z  AI  =  (r  b)a  Y,  Am  , 

1=1  1  ^  1=1  ^ 

which  must  hold  for  all  n.  Therefore, 

AI  =(  rb)'-Am,  . 

'a.  1 


(8) 


(9) 


(10) 
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A  cylindrical  shape  was  chosen  for  the  weights  with 
radius  and  length  and  with  axis  at  the  desired  c.g. 
location.  The  mass  polar  moment  of  Inertia  of  the  ith  weight 
per  unit  span  about  the  rotational  axis  Is  given  by 

Am. 

Equation  11  combined  with  equation  10  shows  that: 


R.  =  b  C2(r2  -  x^)  1^/^.  (12) 

1  a  a  '  ' 

But  the  right-hand  side  Is  Independent  of  1.  Therefore^  the 

subscript  may  be  dropped  from  R.  The  length  of  the  1th  weight 

Is  governed  by  the  required  change  In  massj  as  follows* 


(13) 


where  Is  the  mass  density  of  the  second  set  of  weights. 
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III. 


THEORETICAL  ANALYSIS 


A.  GENERAL 

In  hydrofoil  applications^  an  adequate  definition 
of  flutter  properties  can  be  obtained  from  a  study  of  the 
stability  of  very  small  motions.  Small  motions  should  be 
stable;  otherwise,  they  lead  to  fatigue  conditions  If  not  to 
larger  destructive  motions.  In  practice,  the  large  motions 
usually  are  stable  if  the  small  motions  are  stable.  There¬ 
fore,  time  dependence  can  be  assumed  to  be  proportional  to 
e^^,  where  A.  Is  a  complex  number;  all  other  motions  can  be 
built  up  by  superposition.  (The  procedure  for  superposing 
the  lift  response  .Cor  more  complicated  motions  Is  given  In 
Appendix  A.) 

Luring  flutter,  the  critical  flutter-speed  Is 
associated  with  a  frequency  of  sustained  vibratory  motion. 

The  determination  of  the  critical  flutter-speed  becomes  a 
simpler  problem  v/hen  simple  harmonic  motion  Is  assumed;  It 
Is  easier  to  describe  mathematically  the  hydrodynamic  loads 
for  simple  harmonic  motion  than  It  Is  to  describe  these  loads 
for  more  general  motions.  Therefore,  to  simplify  the  problem, 
the  speed  and  frequency  required  for  sustained  simple  harmonic 
motion  should  be  determined,  rather  than  the  response  of  the 
foil  systjm  as  a  function  of  speed.  In  this  Investigation, 
all  time  dependent  terms  for  simple  harmonic  motion  were 
assumed  proportional  to  e*^^^  (co  real),  rather  than  e^^ . 

Figure  8  shows  the  coordinates  and  nomenclature  used  In  this 
investigation . 

B.  EQUATIONS  OP  MOTION 

The  equations  of  motion  for  the  experimental  model 
were  derived  from  Lagrange’s  equations.  A  lumped-parameter 
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system,  with  negligible  damping  due  to  friction  or  structural 
deformations,  was  assumed  representative  of  the  model. 

Figure  9  shows  the  dynamic  system  and  the  following  defini¬ 
tions  were  used: 

iri^  =  mass  of  translating  parts  per  unit  span; 

m^  =  mass  of  rotating  parts  per  unit  span; 

I  =  mass  polar  moment  of  Inertia  of  rotating  parts 
^  about  Its  center  of  gravity  per  unit  span. 

The  Lagranglan  for  the  translation  and  rotation  displacements 

C  h(  t)  and  a(  t)  1  Is 

I  =  iv"  +  -  IV  -  Iv- 

By  substituting  £  Into  Lagrange's  equation,  the  equations  of 
motion  for  the  system  are 

mh  +  m^x^ba  +  K^h  =  (l5) 

and 

=  «o'  < 

Where  the  dots  Indicate  differentiation  with  respect  to  t  and 

the  following  definitions  have  been  used: 

m  =  m^  +  =  total  oscillating  mass  per  unit  span, 

I  =  I  +  (x  b)^m  =  mass  polar  moment  of  inertia  of 
^  ^  a  a  rotating  parts  about  the  rotational 

axis  per  unit  span. 

and  Q^-,are  the  generalized  force  and  moment,  respectively, 
acting  on  the  system.  By  considering  the  work  done  by  external 
forces  during  a  virtual  displacement  (6h  and  6a),  and 
became  the  hydrodynamic  lift  and  moment  per  unit  span  at  the 
rotational  axis. 

The  equations  of  motion  can  be  put  in  non-dimensional 
form  by  dividing  equation  15  by  TipbU^  and  equation  l6  by 
Trpb^U^  and  by  using  the  following  non-dimensional  parameters: 
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M  =  m/irpb^ 
^  =  m  /m 


=  density  ratio, 

=  coupling  mass  ratio, 


U  ^  =  K.  b^/mU^  =  dimensionless  uncoupled  natural 
^  ^  frequency  In  translation. 


Q,  ^  =  K  b^/l  1]^=  dimensionless  uncoupled  natural 

frequency  In  rotation, 

r  ^  =  I  /mb^  =  dimensionless  radius  of  gyration  about 

rotational  axis. 


s  =  Ut/b 
h(  s) ,  a(  s) 


=  dimensionless  time  variable,  and 

=  dimensionless  translation  and  rotation 
displacements 


so  that 


h(  s)  =  ^h(  t)  and  a(  s)  =  a(  t)  . 


By  Introducing  these  quantities  and  by  nondlmenslonallzlng, 
the  equations  of  motion  become 

p  h"  +/3/Jx^a"  +/rQ|^2h  =  h/np  hU^  ( 1' 


pr  "a"  +  p>iJx  h"  +pr  =  M/^p  b^U^  ( 19, 

where  the  primes  Indicate  dll'I'erentlatlon  with  respect  to  s. 


1.  Flutter  Analysis 


The  description  of  the  mode  shape  associated  with 
the  flutter  of  a  hydrofoil  connected  to  a  hull  Involves  the 
superposition  of  the  normal  vibration  modes  of  the  complete 
structure.  When  this  precise  mode  shape  Is  used  In  flutter 
analyses  It  results  In  a  lengthy  computation  because  there 
are  an  Infinite  number  of  normal  vibration  modes. 


Aeroelastlc  computations  have  shown  that  flutter 
speed  and  frequency  can  be  predicted  accurately  by  the 
Rayleigh- Rltz  procedure.  With  this  procedure,  the  mode  shape 
Is  approximated  by  the  superposition  of  a  finite  number  of 
pre-asslgned  or  assumed  modes.  The  choice  of  modes  Is 
governed  by  the  structure  Involved  and  by  the  results  desired. 
Reference  7  discusses  the  choices  available  and  the  criterion 
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for  their  use.  An  accurate  and  simple  procedure  Is  to  use 
two  or  three  primitive  mode  shapes  that  are  compatible  with 
the  boundary  conditions.  The  assumed  modes  method  of  flutter 
analysis  has  been  compared  with  a  more  exact  treatment 
(reference  3)*  The  two  methods  reasonably  agreed  In  the  range 
of  parameters  used  In  hydrofoil  design. 

In  this  Investigation,  the  assumed  modes  method 
exactly  described  the  flutter  mode  of  vibration;  the  two 
degrees  of  freedom  of  translation  and  rotation  In  the  plane 
of  flow  are  the  only  appreciable  displacements  of  the  model. 

A  linear  combination  of  these  displacements  exactly  des¬ 
cribed  the  flutter  mode-shape.  Thus,  no  approximation  was 
used  In  the  application  of  the  assumed  modes  method. 

In  this  case.  It  was  more  convenient  to  divide  the 
equation  of  motion  by  Trpb^cu^  and  T\ph%^  rather  than  use  equa¬ 
tions  l8  and  19-  Introductlng  the  assumption  of  simple  har¬ 
monic  motion  In  the  form 

h  =  h  be*^^^  and  a  =  a  e*^^^  (20) 

0  0 

and  using  the  quantities  p,  x^,  and  r^,  the  equations  of 
motion  were 

pCl  -  (- — )^(— )^1  h  e"^^  +  p/3x  a  e"^^  =  L/mpb^oo^  (2l) 

and 

PiSx  h  +  pr  2  CP  _  =  MApb'^oo^  (22) 

At  this  point,  an  explicit  representation  of  the 
hydrodynamic  loads  was  Introduced.  The  Theodorsen  results 

cited  In  reference  8  were  used  In  this  analysis.  In  the  design  of 
the  experimental  apparatus,  every  effort  was  made  to  model 
the  assumptions  Involved  In  the  ttieory-- that  Is,  large  end- 
plates  were  provided  with  small  clearances  to  preserve  two- 
dimensional  /'low;  the  .foil  supports  were  designed  so  that 
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the  desired,  "rigid  -  body  motions  were  achieved  at  every 
point  along  the  span  of  the  foil;  no  free  surface,  cavitation, 
ventilation  or  sidewall  effects  were  permitted.  In  addition, 
the  foil  motions  were  assumed  to  be  sinusoidal;  however,  the 
recorded  motions  were  checked  only  visually. 

The  representation  of  the  hydrodynamic  loads  was 
taken  from  reference  9>  where  the  expressions  for  the  hydro- 
dynamic  lift  and  moment  are  given  as 

L  =  TrpbV  +  CL„  -(1  +  a)  (23) 

and 

M  =  Trpb^a)^  C M.  -(^  +  a) L,  1  h  +  CM  -(l  +  M,  )(i  +  a) 

|h  2  h  0  a  a  h  2 

+  (i+  a)  =  :aj  (24) 

In  these  relations,  ab  Is  the  distance  from  midchord  to  the 
rotational  axis  and  Is  positive  If  the  rotational  axis  Is 
aft.  The  coefficients  L^,  L^,  and  are  complex  func¬ 
tions  of  the  reduced  frequency  ( k),  which  Is 

k  =  (25) 

A  tabulation  of  these  coefficients  may  be  found  In  references 
9  and  10. 

These  results  were  Introduced  Into  the  equations  of 
motion  with  the  following  notation: 


hr 

=  L, 

Itr 

L, 

hj 

=  K  I 
hJ 

Lar' 

=  L 
ar 

+ 

a) 

hr 

11 

+ 

a)L,  . 
hj 

M.  ’ 
hr 

=  M, 
hr 

-'I 

+ 

a)  L, 
hr 

Mu 

hj 

=  M,  , 
hJ 

+ 

f— 
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M 

ar 

T  - 

M 

ar 

+  a){L^^ 

+  V 

+ 

+  a) 

^L, 

hr 

»  = 

M 

aj 

+  a)(L^j 

+ 

+  a) 

^L,  , 
hJ 

d  = 

kr 

'M 

aj 

’  + 

1 

L  'M  ’ 

hj  ar 

-  kr' 

“hJ 

f 

kj 

'\r 

e  = 

kr 

'  - 
ar 

L  ’M  ' 

hJ  "aj 

-  kr' 

kr 

‘  + 

'  kj 

f  = 

kj 

'  + 

1  = 

^ar 

»  + 

V 

(26) 

The  subscripts  r  and  J  denote  the  real  and  Imaginary  parts  of 
the  corresponding  complex  numbers.  Equations  23,  24,  and  26 
were  Introduced  Into  equations  21  and  22  and  which  Is 

common  to  all  terms  was  cancelled,  yielding: 

(  CU  (JO  1  (  ) 

L  c  1  -  (^)  ■  h^+  .  0 

(a  )  <  (27) 


and 


CD 


jKPx„+MhP  +  JMhj'[  h„+  a  =  0. 

These  are  the  equations  of  motion  that  the  unloiowns  h  ,  a  , 

^  0  0 

U,  and  (JO  must  satisfy. 


Equations  27  and  28  are  complex,  homogeneous,  simul¬ 
taneous  linear  algebraic  equations  In  the  unknown  amplitudes 
h  and  a  .  A  nontrivial  solution  for  these  unknowns  exists. 

If  and  only  If  the  determinant  of  their  coefficients  vanishes. 
In  this  case,  the  determinant  Is  complex  and  results  In  two 
simultaneous  equations  by  setting  the  real  and  Imaginary  parts 
of  the  expanded  determinant  equal  to  zero  separately.  The 
solution  Is  thus  reduced  to  finding  the  eigenvalues  for  U 
and  03  which  satisfy  two  algebraic  equations. 


However,  U  appears  In  equations  27  and  28  only 
Implicitly  In  the  argument  (k)  of  the  hydrodynamic  coefficients. 
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Because  these  coefficients  are  most  easily  handled  In  tabu¬ 
lar  form  rather  than  functional,  the  argument  usually  Is  con¬ 
sidered  as  an  Independent  parametric  variable.  The  unknown 
frequency  and  one  other  parameter  are  then  considered  the 
dependent  unknowns  and  solutions  to  the  expanded  determinant 
are  found  at  chosen  values  of  k.  The  second  unknown  is 
usually  chosen  so  that  one  of  the  equations  resulting  from 
the  expansion  of  the  determinant  Is  linear  In  both  dependent 
unknowns.  One  method  of  solution  of  the  flutter  determinant 
Is  the  density  variation  method.  Alternate  means  of  solution 
are  described  In  reference  7* 


2.  Density  Variation  Method 


The  density  variation  method  was  particularly 
appropriate  here,  because  the  density  ratio  was  varied  while 
other  parameters  were  kept  constant.  Thus,  the  density  ratio 
and  the  square  of  the  unknown  frequency  ratio  were  taken  as 
the  dependent  unknowns  In  solving  the  flutter  determinant, 
which  Is: 


03,  03^ 

/£!  -  (— 3  «  +  JL.  ' 

03^  03  hr  hj 


03 


^  "hr'  +  J"hj' 


a\  5 


^  +  «ar'+JV 


=0 


(29) 


By  expanding  and  setting  the  Imaginary  part  of  the  resulting 
equation  equal  to  zero,  the  following  linear  relation  exists 
between  the  chosen  dependent  unknowns: 


i=  V{— ) 


03 


^  a  ^hj  '‘03''  d  "aj  a 


(30) 


d  and  f  ar,e  defined  In  equation  26.  By  setting  the  real  part 
of  the  expanded  determinant  equal  to  zero,  and  eliminating 
l/p  by  Introducing  equation  30, 
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(jO  /  0)  0)  1 

(_aw.  2(Jl)a  _  ^  c(J1)2m  t  +  r/L  >]  +  eH^ 

'  oj  ^  j^a  'o)  ''  ar  a  hr  I 

(  tt  a  •  ) 


CO 


to 
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•a  ^hr 


-  «  «ar'  +  '■a’^^hr'  ' 


+  r/  -  (PxJ^  +  J  CM^^'  +  -  fix^ll  +  ej^  =  0, 


(31) 


which  contains  only  one  unknown  (co^/co)^,  in  quadratic  form. 
The  expressions  for  H  and  J  are  defined  by  equation  30  so 
that 


1 


(32) 


At  chosen  values  of  k,  equation  31  is  solved  for  values  of 
(to^/co)^,  which  are  substituted  Into  equation  30  to  determine 
p.  Finally,  the  corresponding  values  of  the-reduced  speed 
are  found  from  the  following  relationship: 


IL  _  l/l^ 


(33) 


This  method  of  flutter  analysis,  Including  appro¬ 
priate  elastic  and  dynamic  assumptions,  Is  believed  to  be 
most  useful  In  hydroi'oll  design  procedures.  The  results  of 
this  analysis  applied  to  the  expei-lmental  conditions  are 
shown  In  Figure  10. 

The  Theodorsen  representation  of  the  hydrodynamic 
forces  limited  this  analysis  strictly  to  simple  harmonic  motions 
In  two  degrees  of  freedom.  If  other  motions  are  Involved, 
which  are  either  not  simple  harmonic  or  other  degrees  of 
freedom,  this  analysis  Is  not  applicable.  Subject  to  these 
restrictions,  the  curve  of  reduced  flutter  speed  versus 
density  ratio  obtained  by  the  density  variation  method  defines 
the  condition-- that  Is,  speed  ( u)  and  frequency  (a))--that 
must  prevail  at  flutter  for  a  particular-  hydrofoil  section. 

At  speeds  above  critical  flutter-speed,  the  motion  of  the 
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foil  will  be  oscillatory  divergent  and  at  lower  speeds  the 
motion  will  be  convergent.  The  curve  plotted  In  Figure  10 
is  the  boundary  between  unstable  motion  at  higher  speeds  and 
stable  motion  at  lower  speeds. 

3,  Stability  Analysis 

In  the  stability  analysis^  the  time  dependence  of 

xt 

the  motions  again  was  assumed  to  be  e  .  Therefore,  from  the 
equations  of  motion,  a  characteristic  equation  was  derived, 
which  was  a  polynomial  In  X.  The  characteristic  values  (the 
roots  of  the  polynomial)  were  found  and  a  quantitative 
measure  of  stability  was  obtained  from  considering  the  mag¬ 
nitude  and  sign  of  the  real  part's  of  all  roots.  In  this 
sense,  the  assumed  modes  method  of  flutter  analysis  Is  not  a 
stability  analysis  because  of  the  simplifying  assumption  of 
sinusoidal  motion. 


In  this  analysis,  the  equations  of  motion  in  non- 


dimensional  form  were  used  : 


/ih"  +  /i/3x^a"  +  -  L/TipbU^  =  0 


/ir  ^a"  +  /i,Sx  h"  +  pr  -  MApb^U^  =  0.  ( I9) 

(X  a  a  a  '  '  \ 


The  hydrodynamic  terms  arp  introduced  by  the  Duhamel  superposi¬ 
tion  Integral  and  the  Wagner  function  0(  s) ,  which  gives  the  time 
dependence  of  the  circulatory  lift  response  to  a  unit  step 
change  In  angle  of  attack.  The  development  of  these  expres¬ 
sions  and  the  following  results  are  described  In  Appendix  A 
of  this  report  and  la  references  7  and  10.  The  Theodorsen 
representation  of  the  lift  and  moment  Is  net  applicable  to 
this  analysis  ( Appendix  b) . 


The  hydrodynamic  force  and  moment  are 
-L/vTPbU^  ='h"  +  a'  -  aa"  +  2W(  s) 


and 

=  ah"  +  +  a)a’  -f  +  a^)n,"  -2(i  +  a)w(s) 


(34) 

(35) 
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In  these  equations,  the  quantity  W( s)  Is  defined  as 


-s 

W(  s)  =H’(o)(/)(3)  +1  0(  s  -  a)H"(a)da, 


(36) 


where 

H»(s)  -  h’  +  a  +  (i  -  a)a'  (37) 

and  where  H'(o)  Is  the  value  of  ]l’(s)  at  s  =  0.  Equation  37 
Is  Identical  to  the  expression  for  the  vertical  downwash 
velocity  at  the  3/^t  chord  point  of  the  I'oll. 


An  approximate  form  of  the  Wagner  function  was 
used  In  this  analysis.  The  following  are  two  approximations 
commonly  used  In  aeronautical  literature  (reference  7) • 


or 


0(s)  ^  1 


-  0.335e"°‘^^ 


0(  s) 


N  s  +  2 
o~+~7" 


(38) 

(39) 


The  former  was  used  because  of  Its  adaptability  to  the 
Laplace  transi'orm  method  of  solution.  Introducing  these 
results  Into  the  equations  ol'  motion,  taking  the  Laplace 
transform,  and  using  the  Initial  conditions  (Figure  8), 


h'(0)  =  a>(0)  =  0,  h(0)  =  h^  and  a(  o)  =  a^,  ( 4o) 

results  In  the  following  pair  of  simultaneous,  linear, 
algebraic,  nonhomogeneous  equations  In.  the  transforms  of  the 
motion  h(p)  and  a(p)  In  which  p  Is  the  complex  Laplace  trans¬ 
form  variable; 

h  |p^(p  +  l  +  20)  +  F'W,  ■-  I H-  a  -  a  -I-  2(4  -  a)0]  + 

[  n  j  I  cx  L. 

p(l  +  20)  I  =  h^p( //, -b  1 T  20)  +a^  jl  +  pCp/^x^ -a-l-2(i-a)03| 

(41) 

and 
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hp^  I  ^  j  ■*■  H  ■'■  ^^  ^ 

+  pC^  -  a  -2(i  +  a)03  +  ^  !  =  h  pCpte  -  a  -2(i  +  a)0] 


+  a 


o  2 


-  a  +  pC/ir 


a  a 

; 

=  -  2( 


+  i  +  a  -  ai-jf 


0"  '  '  a 
1  -2 


-  a*-)^  3  [  , 


(42) 


where  the  bar  Indicates  Laplace  transformation.  The  trans¬ 
form  of  equation  38  yields 

/NTZC/r  r\  --i  rr 

(43) 


0(p)  =  i  - 


p  ”  p  +  0.0455  "  P  +  0.3 


Finally,  by  solving  equations  4l  and  42  for  the  transformed 
responses  h(p)  and  a(p),  then  clearing  the  fraction  yields 


"  Ifpl"  '  ifpf 

where  P,  Q  and  R  are  polynomials  In  p.  ■  and  R  are  of  order 
five  and  Q  Is  of  order  six,  which  permlti-d  the  following 
Inverse  transformation  to  be  used: 


■1  PU 


Q(P 


=  E 


6  P(pJ 


n' 


P  s 
,  n 


(“IS) 


Where  p^  are  of  the  roots  of  the  characteristic  equation 

Q(p)  0  (46) 

and  where  Q'(p  )  Is  the  derivative  of  Q(p)  with  respect  to  p 
n 

evaluated  at  p  =  p^.  These  results  may  be  confirmed  In  any 
textbook  on  operational  methods  (reference  ll) .  An  Iterative 
process  based  on  Newton's  method  was  used  to  find  the  quadratic 
factors  of  equation  46  and  Is  described  In  reference  12.  The 
Interest  here  lies  In  the  complex  roots  of  equation  46  and  In 
particular  the  root  In  which  the  real  part  becomes  positive 
above  some  critical  speed.  The  results  were  obtained  In  the 
form  of  the  overall  damping  ratio  ( a/co^) ,  which  Is  the  ratio 
of  the  real  part  of-  kto  the  uncoupled  rotational  frequency. 

If  the  complex  root  of  equation  46  Is 

p  =  u  +  1  V  (47) 
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then  the  overall  damping  ratio  Is  found  from  the  real  part  of 
p  and  the  frequency  ratio  f^om  the  Imaginary  part  as 

follows 


The  overall  damping  ratio  as  predicted  by  this  method  for  the 
experimental  conditions  Is  plotted  In  Figures  11  through  I7 
with  the  corresponding  experimental  results  and  the  frequencies 
are  presented  In  Figures  18  through  25. 


C.  CRITICAL  DENSITY  RATIO 


An  examination  of  the  theoretical  results,  (parti¬ 
cularly  those  shown  In  Figure  lO)  will  show  that  the  reduced 
flutter  speed  becomes  Infinite  at  some  value  of  density  ratio. 
This  value  can  be  predicted  by  using  the  following  Information: 
u/bo)^  03^  when  k-»0  for  finite  values  of  An  exten¬ 

sion  of  the  assumed  modes  analysis  leads  to  an  explicit  equa¬ 
tion  for  Up/boD^  when  these  limits  are  introduced.  This  result 
Is  developed  In  reference  3>  where  the  equation  Is  shown  to  be 
a  ratio  of  two  polynomials.  Setting  the  denominator  of  this 
result  equal  to  zero  yields  the  relation  between  the  para¬ 
meters  that  must  be  satisfied  at  the  asymptote ■ where 
Up/bcc^  00  and  m  the  critical  density  ratio.  It  Is 

found  that  /i_„  decreases  with  co,  /u)  and  a  very  simple  rela- 
On  h  CL 

tlon  Is  found  for  the  limiting  case  =  Oj 


where  x^  +  a  is  the  c.g.  location  measured  in  half-chord 
lengths  from  midchord,  positive  If  the  c.g.  Is  aft.  This 
result  Is  plotted  In  Figure  25- 


By  this  procedure,  an  asymptotic  behavior  Is  pre¬ 
dicted  for  the  flutter  speed  as  p  approaches  However, 

OK 

this  asymtote  Is  approached  from  the  low  p  side  In  many  cases 
so  that  Is  not  always  the  minimum  density  ratio  at  which 
flutter  is  predicted. 
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IV.  DISCUSSION 


A.  GEInIERAL 

In  these  experiments,  flutter  was  induced  under 
controlled  and  somewhat  artificial  conditions  to  validate  the 
application  of  airfoil  theory  to  hydrofoil  phenomena.  Con¬ 
sequently,  these  theoretical  studies  were  carried  out  in  ac¬ 
cordance  with  aeronautical  practices. 

B.  EXPERIMENTAL  PROCEDURE 

At  each  mass  condition  ( ju  ,  where  n  =  0,  1,  2,  4, 

7,  9)  ll)  two  bits  of  experimental  Information  were  obtained; 
the  overall  damping  ratio  ( o/co^) ,  and  the  frequency  ratio 
(oj/cO^)  .  The  model  was  locked  in  position  and  when  the  apparatus 
attained  the  desired  speed,  the  model  was  released  and  the 
resulting  motion  was  recorded.  The  speed  of  each  successive 
test  was  Increased  until  the  motion  of  the  model  was  diverg¬ 
ent,  in  which  case  the  apparatus  was  quickly  brought  to  rest 
to  avoid  destructive  or  violent  motions.  Flutter  was  obtained 
in  each  test  series  except  series  0  ( =  0,758). 

A  time  signal  was  recorded  simultaneously  with  the 
motions  to  measure  frequency  (cu)  .  Because  the  model  had  two 
degrees  of  freedom,  two  values  of  co  were  obtainable ,  The  coupled 
zero-speed,  rotation  mode-shape  led  to  flutter.  The  dairiplng 
in  the  translation  mode  Increased  rapidly  with  speed,  con¬ 
veniently  aiding  Its  elimination  from  the  records.  By  properly 
choosing  Initial  conditions,  appreciable  response  in  the 
translation  mode  was  eliminated. 

The  second  piece  of  Information,  the  overall  damping 
ratio  (a/cu^),  was  obtained  using 
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(49) 


a 

cr  1  /  OD  X  p  -tfi/  1  +  1\  -1 

CD  "  aT'cD  1  ^  a  ^  ’ 

a  a  0^ 

where  1  Is  the  number  of  cycles  analyzed.  The  quantity  In 
square  brackets  Is  proportional  to  the  slope  of  a  straight 
line  faired  through  the  successive  amplitudes  of  oscillation 
plotted  on  semi-log  graph  paper. 

The  tests  were  performed  In  two  groups  ( series  k, 

7j  11  and  series  0,  1,  2,  9,  respectively) ,  separated  by  two 
weeks,  during  which  time,  some  of  the  data  were  analyzed. 
During  this  time  It  was  r-eallzed  that  the  rotation  mode  of 
vibration  was  leading  to  flutter;  subsequently,  the  model  was 
held  at  an  Initial  angle  of  attack  rather  than  with  an  Initial 
translation  displacement,  as  In  the  I'lrst  group  of  tests. 
During  the  second  group  of  tests,  an  Investigation  of  the 
effect  of  Initial  conditions  on  the  experimental  measuiements 
showed  that  the  results  of  the  i'lrst  group  could  still  be 
used,  because  of  the  large  amount  of  damping  In  the  transla¬ 
tion  mode  of  vibration;  by  repeated  runs  at  a  particular 
speed,  the  results  were  found  to  be  reproducible.  (See  points 
at  U/bCD^  =  1.11  In  Figure  14  oi'  at  U/bCD^  =  1,21  In  Figure  17.) 

C.  EXPERIMENTAL  RESULTS 

The  maximum  speed  obtainable  was  considerably  lower 
than  anticipated.  The  apparatus  was  designed  on  the  basis  of 
a  maximum  speed  of  30  i’eet  per  second.  Actually,  this  speed 
was  15  feet  per  second  ( u/bco^  =  I.67,  for  series  0,  p  =  O.758, 
where  the  limitation  was  most  severe  and  the  only  mass  condi¬ 
tion  In  which  the  maximum  speed  was  obtained).  The  experimental 
results  i'or  the  present  Investigation  are  somewhat  Inconclu¬ 
sive  In  this  regard:  It  Is  not  certain  that  flutter  could  not 
be  Induced  at  a  higher  speed  Cor  series  0  (Figure  ll) .  How¬ 
ever,  a  sharp  upward  turn  In  reduced  Cluttei’  speed  does  occur 
as  p  decreases,  (Figure  lo) ,  whlcli  at  least  partially  sub¬ 
stantiates  the  exlstance  of  an  asymptote  as  predicted. 
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D 


THEORETICAL  RESULTS 


The  analysis  by  the  two  theoretical  methods  that 
can  be  used  to  predict  speed  and  frequency  at  flutter  has 
shown  that  ( l)  the  assumed  modes  method  of  flutter  analysis 
leads  to  a  prediction  of  at  which  the  flutter  speed 
Increases  asymptotically  and  (2)  the  superposition  method  of 
stability  analysis  provides  Information  for  speeds  other  than 
the  flutter  speed;  these  results  were  obtained  In  addition  to 
the  prediction  of  critical  flutter- speed  and  frequency. 

The  flutter  speeds  as  predicted  by  the  assumed 
modes  analysis^  which  Involved  an  exact  modal  description  In 
this  casOj  are  shown  In  Figure  10  and  the  corresponding  fre¬ 
quencies  for  each  test  condition  are  shown  In  Table  II.  The 
critical-density  ratio  as  predicted  by  this  analysis  at 
oOp^/o)^  =  0  Is  plotted  In  Figure  25*  The  difference  between 
series  0^  1,  2,  9  and  series  1,  11  (Figure  lO)  is  caused 

by  the  change  in  x  between  these  two  groups^  thus  causing  a 
I  ^ 

slight  shift  In 

On 

Figures  11  through  I7  show  the  overall  damping  ratio 
obtained  from  the  superposition  analysis,  which  reduces  to  the 
well-known  critical  damping  ratio  for  a  one-degree  of  freedom 
system  vibrating  freely  In  air,  plotted  against  the  reduced 
speed  for  each  mass  condition.  The  predicted  flutter  speed  is 
found  where  u/co^  goes  to  zero.  A  comparison  of  these  speeds 
with  the  theoretical  results  from  the  assumed  modes  analysis 
In  Figure  10  for  corresponding  p's  shows  the  anticipated 
Identical  agreement  between  the  two  methods.  The  use  of  the 
approximation  In  equation  38j  thus  agrees  well  with  more  exact 
theory.  No  flutter  Is  predicted  for  series  0  or  1.  Figures  I8 
through  2^1  show  the  frequency  ratio  versus  reduced  speed,  as 
predicted  by  the  supei'posltlon  method. 

Other  results  of  the  superposition  analysis  are 
tabulated  In  Table  III  for  series  9  (p  3.0-)). 


The  theoretical  results  described  here  were  obtained 
on  an  electronic  digital  computer.  The  program  for  both  the 
assumed  modes  and  the  superposition  analysis  have  been 
retained  for  future  use. 
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V. 


THEORY  COMPAIffiD  WITH  EXPERIMENT 


The  results  of  series  9j  and  11  show  relatively 
good  agreement  between  theory  and  experiment  for  the  predicted 
flutter  conditions  (Table  II) .  As  cited  in  reference  4  and 
as  indicated  in  Figure  10.,  the  theory  gave  a  conservative 
estimate  of  flutter  speed  for  these  conditions.  No  apparent 
reason  could  be  found  to  doubt  the  results  of  series  7 
(p  =  2.o6)  in  Figure  10;  however,  there  is  nothing  to  cause 
these  results  to  be  out  of  line  with  the  other  series.  A 
lower  experimental  value  for  the  reduced  flutter  speed  at 
this  mass  condition  would  be  in  agreement  v^ith  the  antici¬ 
pated  experimental  results. 

Between  series  7  and  'i,  the  Inadequacy  of  the  theory 
to  predict  the  flutter  speed  with  accuracy  becomes  more 
apparent  and  the  resulting  error  becomes  more  dangerous. 
Although  the  good  agreement  in  frequency  ratio  is  retained  at 
the  lower  density  ratios,  the  overall  damping  ratio  is  con¬ 
siderably  overestimated- b-y  the  theory;  in  fact,  at  series  1 
the  theory  predicts  freedom  from  flutter  where  flutter  was 
obtained  during  the  experiments.  In  the  case  of  series  0, 
flutter,  though  It  did  not  occur,  was  Imminent;  the  theory 
predicts  no  flutter  at  all  speeds.  Even  in  this  circumstance 
the  frequency  ratio  is  well  predicted  (Figure  l8)  and  is 
within  9  percent. 

The  results  of  reference  ^1,  indicated  a  discrepancy 
between  theory  and  experiment  for  the  prediction  of  flutter 
speed  at  low  density  ratios.  However,  tliey  were  not  able  to 
extend  their  investigations  to  a  sufficiently  low  value  of 
density  ratio  to  check  the  existence  of  an  asymptote  as  pre¬ 
dicted  by  equation  ^19.  Tliese  result. s  are  sununarlzed  in 
Figure  1.  At  high  values  of  density  ratio  (p  --  3-0)  tlie 
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theoretical  and  experimental  results  were  essentially  in 
agreement j  or  the  theory  was  conservative. 


The  present  results  shown  in  Figure  10  substantiate 
the  results  of  reference  ^1 .  In  addition,  the  experimental 
results  at  series  0( p  =  0.758)  strongly  Indicate  the  exist¬ 
ence  of  an  asymptote  In  flutter  speed.  There  Is  definitely  a 
sharp  upward  trend  shown  as  p  decreases.  The  position  of 
this  asymptote  Is  not  correctly  predicted  by  theory.  The 
theoretical  and  experimental  value  of  p^^^  are  plotted  In 
Figure  25-  The  comparison  Is  somewhat  worse  than  Is  shown 
there  since  the  theoretical  curve  Is  derived  for  the  case 


=  0  and  In  reference  3^ 

with  01  /oj  . 

”r  a 


It  is  shown  that  Increases 


Some  more  definite  InI'ormatlon  about  the  discrepancy 
can  be  gotten  from  the  results  of  the  stability  analysis. 

The  experimental  and  theoretical  overall  damping  ratios  are 
compared  in  Figures  11  through  17.  There,  It  Is  seen  that 
a/o)^  is  overestimated  by  the  theory  up  to  series  7  (p  =  2.08) 
and  from  that  point  on  the  predicted  overall  damping  ratio  is 
less  than  that  obtained  In  the  experiments.  The  frequency 
ratio  Is  very  well  predicted  throughout  as  shown  In  Figures  l8 
through  24.  We  have  found  that  the  frequency  ratio  depends 
mainly  on  terms  proportional  ro  h(  t)  and  a(t)  or  their  second 
derivatives;  however,  the  damping  ratio  depends  mostly  on 
terms  proportional  to  the  first  derivatives.  The  latter 
contain  the  circulatory  responses-- that  is,  the  Theodorsen 
or  Wagner  functions;  tlie  former  are  not  as  sensitive  to  these 
functions.  Since  the  fi'equency  ratio  Is  fairly  well  predicted 
and  the  damping  Is  not.  It  may  be  surmised  tliat  the  effective 
mass  and  Inertia  terms  in  tlie  theory  are  more  accurate  than 
the  circulatory  terms.  Othci'  Investigators  have  used  this 
hypothesis  as  mentioned  in  reference  j. 


VI. 


CONCLUSIONS 


1.  A  conservative  prediction  of  flutter  speed  may  be 
expected  at  hlgh-denslty  ratios. 

2.  The  theory  does  not  correctly  pi'edi.ct  the  flutter 
speed  near 

3.  The  theory  does  not  correctly  predict  the  location 
of  Mqpj- 

The  frequency  Is  very  v^ell  predicted  at  flutter  and 
at  lower  speeds  at  all  density  ratios. 

5.  The  overall  damping  ratio  does  not  agree  with 
experimental  results. 

6.  The  circulatory  terms  In  the  theory  are  more  liable 
to  doubt  than  the  effective  mass  or  spring  terms. 


VII. 


RECOMMENDATIONS 


Ii\arther  investigations  are  In  order  to  describe 
more  fully  the  discrepancy  that  has  been  shown  to  exist. 
Several  different  values  of  c.g.  location  and  radius  of 
gyration  should  be  used  In  further  Investigations.  In  addi¬ 
tion  to  the  results  presented  here,  an  experimental  determi¬ 
nation  of  amplitude  ratio  and  phase  angle  should  be  attempted. 
Single  degree  of  freedom  studies  In  water  should  be  carried 
out,  to  provide  separate  comparisons  between  theory  and  exper¬ 
iment  for  each  component  of  the  lift  and  moment  expressions. 
Also,  zero-speed  responses  must  be  analyzed  to  check  the  pre¬ 
diction  of  added  mass,  added  Inertia,  and  added  mass  coupling 
terms . 
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TABLE  I 


PROPERTIES  OF  MODEL  TENITE  II,  FORMULA  233.  FLOW  MS 


Flow  Temperature  284°F 

Specific  Gravity  1.22 

Modulus  of  Elasticity  1.30  x  lO^psl 

Deformation  under  Load  (Cold  Plow)  33% 

Water  Absorption  (24  Mrs.  Immersion) 

Total  Weight  Gained  1.3% 

Soluble  Matter  Lost  0.1% 


Accelerated  Aging  Weight  Lost 

(72  hrs  at  i8o°F)  1.6% 


TABLE  II 

PARAMETERS  AND  FLUTTER  CONDITIONS 


m  X  10 


Series 

No. 

Ji 

X 

a 

2 

r 

a 

(-^) 

(j) 

a 

U  /bo.' 

F  a 

exp.  the . 

0) 

F 

exp. 

the . 

IK  2 

lb  sec 

in^ 

a 

rad 

sec 

0 

.758 

.496 

.900 

.251 

-- 

-- 

-- 

2.00 

36.0 

1 

.883 

.496 

.895 

,264 

1.22 

-- 

,960 

-- 

2.33 

32.3 

2 

1.008 

.496 

.891 

,260 

1.14 

1.84 

.970 

1 ,000 

2.66 

30.2 

4 

1.285 

.  465 

.862 

,269 

1.08 

1.42 

.980 

.977 

3.39 

26.8 

7 

2.06 

.465 

.839 

.277 

1.53 

1  ,29 

.910 

1.002 

5.43 

21.3 

9 

3.03 

,496 

c 

,258 

1  .44 

1  .32 

.980 

1  .040 

8.00 

17.5 

11 

4.07 

.  465 

.896 

.2f.8 

1.55 

1  .45 

.^85 

1.020 

10.74 

15.0 

b  =■  3  inchns,  =  -0.5 


I 
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REDUCED  FLUTTER-SPEED 


REGION  I  a)  NO  FLUTTER 

b)  NO  FLUTTER  THEORETICALLY,  NO  EXPERIMENTAL  INFORMATION 
REGION  2  FLUTTER  FOUND  EXPERIMENTALLY,  BUT  NONE  PREDICTED 

REGION  3  THEORY  PREDICTS  FLUTTER  AT  HIGHER  SPEED  THAN  THAT  FOUND  IN  EXPERIMENTS 
REGION  4  THEORY  GIVES  CONSERVATIVE  PREDICTION  OF  FLUTTER  SPEED 


DENSITY  RATIO 


FIGURE  I.  PREVIOUS  EXPERIMENTAL  RESULTS 


FIGURE  4.  APPARATUS,  PORT  SIDE 


I 


EQUILIBRIUM 


h(0)»  hpb  ,  a(0) »  Qq 

h(0)*  a(0)“  0 


FIGURE  8,  REPRESENTATIVE  HYDROFOIL  ORIENTATION  FOR 
THEORETICAL  ANALYSIS 


FIGURE  9.  SCHEMATIC  DIAGRAM  OF  SYSTEM 
USED  FOR  DYNAMIC  ANALYSIS 


SERIES  O 


J _ I _ I _ \ _ L 

IT)  ^  lO  rj  - 

°(nq 


a33ds  a3iimj  a3Dna3a 


FIGURE  10.  REDUCED  FLUTTER- SPEED  VERSUS  DENSITY  RATIO 


ERSUS  REDUCED  SPEED  (STABILITY  ANALY  SIS,  SERIES  O./i  =  0.758) 


I> 


I 
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FIGURE  12.  OVERALL  DAMPING  RATIO  VERSUS  REDUCED  SPEED  (STABILITY  ANALYSIS,  SERIES  \,p.  =  0.883) 


ILITY  ANALYSIS,  SERIES  2,  =  1.008) 


FIGURE  14.  OVERALL  DAMPING  COEFFICIENT  VERSUS  REDUCED  SPEED  (STABILITY  ANALYSIS  SERIES  4 
1.285)  : 


FIGURE  15.  OVERALL  DAMPING  COEFFICIENT  VERSUS  REDUCED  SPEED  (STABILITY  ANALYSIS.  SERIES  7  u. 

2-08)  ’  ^ 


ED  SPEED  (STABILITY  AN ALYSIS ,  SERIES 


.2  0 


<9) 
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cP 
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FIGURE  18.  FREQUENCY  RATIO  VERSUS  REDUCED  SPEED  (STABILITY  AN  A  LYSIS,  SERIES  O,  /X  =  0.758) 


FIGURE  2a  FREQUENCY  RATIO  VERSUS  REDUCED  SPEED  (STABI LIT Y  ANALYSI S,  SERIES  2,  /i.  =  1.008) 


0  VERSUS  REDUCED  SPEED  (ST> 


FIGURE  23.  FREQUENCY  RATIO  VERSUS  REDUCED  SPEED  (STABILITY  ANALYSIS,  SERIES  9,  /i  =  3.03) 
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APPENDIX  A:  STABILITY  ANALYSIS  BY  SUPERPOSITION 


I.  INTRODUCTION 

To  predict  the  stability  of  small  elastic-deforma¬ 
tions  of  a  hydrofoil  section,  three  types  of  forces  must  be 
related  to  the  motions  of  the  foil:  Inertial  reactions, 
elastic  restraints,  and  hydrodynamic  loads.  The  Inertial 
forces  can  be  described  by  either  Lagrange’s  or  Newton’s 
equations;  to  describe  the  elastic  forces,  strain-displacement 
and  stress-strain  relationships  must  be  established.  (Hooke’s 
law  can  be  used  to  relate  the  stresses  to  the  strains.) 
Hydrodynamic  loads  due  to  sinusoidal  motions  have  been  de¬ 
scribed  mathematically  by  Theodorsen  (reference  8). 

A.  HYDRODYNAMIC  LOADS  DUE  TO  SINUSOIDAL  MOTIONS 
The  sinusoidal  motions  are  expressed  as 

h(s)  =  h^e-^^® 

a(s)  =  (A-1) 

where,  h^  and  are  small  dimensionless  amplitudes  of  heave 
and  pitch,  and  may  be  complex;  k  Is  reduced  frequency;  and  s 
Is  dimensionless  time.  The  two-dimensional  lift  (l)  and 
moment  (m)  per  unit  span  for  a  sinusoidal  motion  are  expressed 
as 

L/7TpbU^  =  -Ch"  +  a’  -  aa"  ]  -  2C(k)lIh’  +  a  +  (i  -  a)a’l 
M/7rpb^U^  =  Cah"  -  (^  +  a-)a’  -  (^  -h  a^)a"]  +  2(^  +  a)c(k) 

Ch’  +  a  -t  (i  -  a) a’  3  (  A-2) 


A-1 


27r 

nj  n  = 


(a-5) 


Let  S  =  TU/b  be  the  dimensionless  period  of  the  motion.  Thus, 
the  Fourier  coefficients  are 


=  -c  I 

o  SJ 


r 


rS 


h( s) dSj 


1 

0  S 


ai 


(3)ds, 


h 


n  S 


-Jk  3 


h(  s)  e  ^  ds,  a  =  ^ 

n  S 


a(  s)  e 


-Jk„3 


n 


ds . 


(a-6) 


To  find  the  lift  response  to  this  motion,  the  response  for  a 
single  Fourier  component  must  be  obtained,  which  Is  expressed 
as 


L^TTphU^  =  -e 


n 


-knX  +  “n  +  + 


{l+  (i-  a)jk^}a„]: 
and  the  moment  response  Is  expressed  as 


(A-7) 


Jk  s  (  , 

M^TTphU^  =  -e  ^  ^n  " 

2(i  +  a)c(k^)Cjk^h^  +  |l  +  (1  _  a)jk^|a^l|  (a-8) 

The  lift  and  moment  response  to  h^  are  zero;  however,  there 
Is  a  steady  response  to  the  mean  angle  of  attack  (cc^): 


L  ApbU^  =  -2a  =  L  ' 

o'  ^  0  0 

M  Apb^U^  =  +2{i  +  a)a  =  M  '  (A-9) 

o'  2  '  0  0  ' 


The  sum  of  the  responses  to  the  components  of  the 
motion  (equation  A-4)  results  In  the  total  lift  and  moment: 


A-3 


and  moment  acting  on  the  foil.  Thus^,  the  case  of  a  foil 
operating  In  oblique  seas  may  be  analyzed.  These  results 
were  derived  assuming  a  rigid  chord-section.  Elastic  defor¬ 
mations  In  camber  could  be  Included  using  the  results  of 
Spielberg*.  The  affects  of  aspect  ratios  as  low  as  2  can  be 
Included  by  using  the  table  of  aerodynamic  coefficients  by 
Relssner  and  Stevens**. 

D.  ARBITRARY  MOTIONS 

The  Theodorsen  results  can  be  extended  to  arbitrary 
motions***.  For  this  motion,  the  nonclrculatory  terms  In 
equation  A-10  are  not  changed  and  only  terms  containing  C!(k^) 
are  considered.  Also,  the  circulatory  lift  and  moment  terms 
depend  on  the  motion  In  exactly  the  same  way-- that  Is,  through 
the  vertical  velocity  of  the  3/^^-chord  point,  which  Is 

H'  (  s)  =  h’ (  s)  +  a(  s)  +  (^  -  a)  a’  (  s)  .  (  A-l4) 

Usjjig  this  observation,  the  circulatory  lift  response  to  a 
single  Fourier  component  of  reduced  frequency  k  Is,  for  a 
unit  amplitude  of  H'(s), 

AL  /TrpU^b  =  -  2C(k)e'^^®  (A-15) 

0 

The  Fourier* Integral  representation  of  an  arbitrary  motion 
H’(s)  becomes 


*Splelberg,  Irvin,  N.,  "The  Two-Dimensional  Incompressible 
Aerodynamic  Coefficients  for  Oscillatory  Changes  In  Airfoil 
Camber, "  Wright  Air  Development  Center,  Technical  Note 
WCNS  52-7,  18  August  1952. 

**Relssn0r,  Eric  and  Stevens,  John  E.,  "Effect  of  Finite  Span 
on  the  Airload  Distributions  for  Oscillating  Wings. 
II-Methods  of  Calculations  and  Examples  of  Application", 

NACA  TN  1195,  19^(7. 

***Wagner,  H.,  "Uber  die  Entstehung  des  Dynamlschen  Auftrlebes 
von  Tragflugeln,  "  Z.  Ange.  W.  Math.  Mech.,  Bd.5,  Heft  1, 
February,  1925* 


Consider,  a  step  change  In  H'(s): 


0  , 


h’(0)  +  a(0)  +  (^  -  a)a’(o),  for  s  >  0 


H'(s)  = 

The  Fourier  Integral  for  such  a  motion  Is 


for  s  <  0; 


H’(s)  =  i  Ch'(O)  +  a(0)  +  (|  -  a)a'(o) 


(A-21) 


dk  (A-22) 


Comparing  equation  A-22  and  equation  A-l6  shows  that 

7i(k)  =  Ch»(0)  +  a(0)  +  (|  -  a)a'(o)  ]  ^  ,  (A-23) 

for  this  motion.  The  circulatory  lift  response  to  this  step 
change  Is  found  by  substituting  equation  A-23  Into  equa¬ 
tion  A-l8: 


L^pU^b  =  -  Ch’(O)  +  a(o)  +  (|  -  a)a’(o)  ]  /  e^^®  dk 

J  *  CO 

(A-24) 

The  Wagner  function  C^(s)  ]  Is  defined  as  the  time  dependence 
of  the  circulator  lift  response  to  a  unit  step  change  In  H*(s) 
at  s  =  0  or 


/oo 

^  eJ^®  dk  (A-25) 

-oo 

which  when  substituted  Into”  equation  A-24  yields 

L  /TTpU^b  =  -  2Ch'(0)  +  a(o)  +  (i  -  a)  a' ( O)  Jj2((  s)  (A-26) 
c  ^ 

Dividing  the  Theodorsen  function  into  Its  real  and  Imaginary 
parts  CF(k)  and  G(  k)  3,  respectively,  results  In  two  simpler 
relations  for  the  Wagner  function: 


(A-31) 


M  =  TTpb^u^  Call"  +  (i  +  a)a»  +  (i  +  a^)  a"  D  + 

^  o 

2'nph^JJ^  (|+  a)  H'(o)^(3)  +j  0(s  -  y)  ^H»(7)d7| 


where 


j2f(s)  = 


F(  k) 


sin  ks  dk  , 


(A- 32) 


0 


,  s  <0; 


H'(s)  =(h'(0)  +  Ua(0)  +  b  (|  -  a)a'(o)  =  H'(o),  s  =  0; 


h'(s)  +  Ua(  s)  +  b  -  a)a'(s) 


,  s  >  0.(A-33) 


II.  ANALYSIS  OF  SYSTEMS  WITH  TWO  LEGREES-OF-FREEDOM 


The  lift  and  moment  equations  were  used  in  the 
stability  analysis  to  determine  the  overall  damping  and  fre¬ 
quency  ratio  associated  with  the  experimental  conditions.  To 
do  thlSj  equations  A-30  through  A-33  are  substituted  into  the 
equations  of  motion,  equations  I8  and  19 .  The  equations  can 
be  written  as 


(iu+  l)h"  +  +  (/ij3x^  -  a)  a"  +  a'  +  2  CH»(o)j?^(s)  + 


7)H"(7)d7:  =  0 


(A-34) 


(/i/3x  -  a)h"  +  (pr^  +  ^  +  a^)a"  +  (i  -  a)a’  +  pr^  + 

Ow  Ob  O  ^  Oj  cc 

2  CH»(0)i?^(s)  +j  ^{s  -  7)H"(7)d7:  =  0.  (A-35) 


A- 9 


Q  _  (1  +  a)fjn  ^ 


{ A-36  cont'd) 


hoSa 


=  h„(M  + 


o  a 


2n  2 


a  iJ-r 
0  ot  <x 


ho«f 

h^Cl  -  (n  +  1)(|  +  a)  -  (pf)x^  -  a)  :  -  -  a)  + 

+  1'^ 


=  h  +  a  [(i  -  -  13 


0  cc  a 


a  Q 

o^a 


(A- 37) 


“0% 


9-) 


'^o^d 


%^f 


Cho  +  (2  -  a)a^D  Q^ 


(A-38) 
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The  polynomials  In  equation  44  are  now  defined  so  that  the 
Inverse  Laplace  transform  Is  Introduced  to  determine  h(  s) 
and  a( s) .  Referring  to  any  text  on  Laplace  transformation 
the  following  theorem  applies*.  If  Q(p)  Is  a  polynomial  of 
degree  1,  with  1  distinct  zeros  p  =  p^,  Pg,...,  p^^^  and  P(p) 
Is  a  polynomial  of  degree  1  -  1  or  less,  then 


0  '  (  Dv, 


(P  -  Pn) 


(A- 43) 


To  each  term  of  equation  A-43  the  following  Inverse  transform 
applies 


1 

P  -  P. 


( A- 44) 


Summing  this  result  for  each  of  the  six  zeros  of  Q(p)  yields 


r-1  W 


6 

=  E 


P(Pn) 


(A- 45) 


The  six  roots  of  the  denominator  lIGi(p)  D  must  be 

determined.  A  process  based  on  Newton’s  Iterative  method  was 
** 

used.  In  this  process,  the  quadratic  factors  of  equation  46 
are  determined  first,  then  the  pairs  of  roots  for  each  quad¬ 
ratic  factor  are  found.  The  polynomial  to  be  factored  Is  of 
sixth  order  so  that  three  quadratic  factors  and  six  roots 
must  be  determined. 

Each  root  Is  complex: 

p  =  u  +  jv 
^n  n  n 


n  =  1,  . . . ,  6 


(A- 46) 


*Hlldebrand,  F.  B.,  ADVANCED  CALCULUS  FOR  ENGINEERS,  Prentice- 
Hall,  Inc.,  New  York,  1954. 

**Descrlbed  In  Scarborough,  J.  B.,  NUMERICAL  MATHEMATICAL 
ANALYSIS,  The  Johns  Hopkins  Press,  Baltimore,  1950. 
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APPENDIX  B:  GENERALIZATION  OF  THEODORSEN'S  FUNCTION 
FOR  CONVERGENT  (STABLE)  OSCILLATIONS 


by  Paul  Rltger 

The  Theodorsen  function  (reference  B-l)was  developed 
for  the  case  of  a  real  argument.  When  this  function  was 
extended  to  Include  complex  arguments,  two  different  and 
Incompatible  conclusions  were  made,  (references  B-2  and  B-3) 
To  resolve  this  discrepancy,  the  following  rigorous  analysis 
was  developed. 

Theodorsen  (reference  B-l)  showed  that  the  clrcula- 
tlonal  part  of  lift  forces  caused  by  translational  motions 
with  unit-magnitude  velocities-- that  Is,  h(  s)  =  e'^^®--ls 

L(s)  =  2p7rUbC(  Oh'(s)  (B-1) 

where  the  Theodorsen  function  C( C)  Is 


The  function  C(  C)  Is  properly  defined  by  equation  B-2,  If  and 
only  If  the  complex  number  C  Is  restricted  so  that 

Iin(d  <  0  (B-3) 


This  Inequality  Is  satisfied  In  the  case  of  divergent  oscil¬ 
lations.  If  Im(  0  =  0,  the  Integrands  In  equation  B-2  behave 
like  sin  ^x  at  «>;  hence,  the  Integrals  would  be  oscillatory 
divergent.  If  Im(  0  >  0,  which  Is  true  for  convergent  oscil¬ 
lations,  then  both  Integrals  In  equation  B-2  become  Infinite. 
Hence,  the  Theodorsen  function,  as  defined  by  equation  B-2, 

Is  meaningless  for  Im(  0  >  0. 


B- 1 


points  then  It  Is  essential  to  specify  which  branch  of  the 
function  Is  to  be  used.  In  the  present  case^the  physically  cor¬ 
rect  branch  Is  partly  determined  by  equation  B-2.  That  Is, 
we  must  be  sure  to  choose  a  branch  of  C^( C)  that  coincides 
with  g(  0  Irn(  d  <  0^  l.e.,  -  tt  <  arg  C  <0.  If  the 
so-called  "principal  branch"  of  K^(jd  Is  chosen,  this  con¬ 
dition  Is  satisfied  and,  hence,  this  seems  to  be  the  most 
natural  choice.  Luke  and  Dengler  have  used  another  branch  as 
will  be  shown  below. 


The  "principal  branch"  of  Is  usually  defined 

(see  Watson,  reference  B-5,  p.  7?)  In  terms  of  the  function 
I  .  That  Is,  In  general,  we  have  the  following  definition  of 
the  many-valued  function  z) 


sin  V7T 


(z  /  o) .  (b-6) 


Now,  I^(z)  Is  In  turn  defined  In  terms  of  J^(z)  and  the 
principal  branch  of  I^( z)  Is  given  by 


I^(z) 


e-Wj/2 

e3v7rj/2 


J^(ze 

J^(ze 


7rj/2)  ^ 
-3'n-j/2 


-7r  <  arg  z  <  7t/2^ 

7r/2  <  arg  z  <  tt.  (B-7) 


In  equation  B-7,  to  be  precise,  the  principal  branch  of  is 
to  be  taken,  which  Is  the  reason  for  a  split  definition 
of  I  (z).  If  z  in  equation  B-7j  lies  In  the  second  quadrant, 
l.e.,  7r/2  <  arg  z  <  tt,  then  ze^"^^  (l.e.  Jz)  would  be  such 
that  77  <  arg  (ze^'^'^^)  <  3'7t/2  but  this  would  then  give  a  value 
of  which  Is  not  on  the  principal  branch.  The  principal 
branch  of  (see  Watson  p.  ^I^l)  Is  obtained  by  restricting 

arg  ^  to  -77  <  arg  C  < 

To  carry  this  all  the  way,  is  defined  by 

00 

j„(  c  =  ( ch  D  a„(  ch 

K=o 


b-3 


(B-12) 


This  definition  of  K  (z)  Is  often  written 

n 

K^(z)  =  J 

which  Is  correct  for  the  multi-valued  function  K^(z)j  but  If 
one  wishes  to  restrict  oneself  to  the  principal  branch,  then 
equation  B-11  must  be  used. 

The  results  of  Luke  and  Dengler  are  obtained  by 
using  equation  B-12  Instead  of  equation  B-11  and  hence  must 
be  Interpreted  with  this  fact  in  mind.  In  particular,  they 
are  concerned  with  the  values  of  C( C)  In  the  vicinity  of 
arg  C  =0.  Hence,  by  equation  B-5  they  use  K^(jC)  near 
v^here  arg  =  0,  l.e.  where  arg  JC  =  1^/2.  So,  if  arg  C 

Is  slightly  greater  than  zero,  they  are  using  values  which 
are  no  longer  on  the  principal  branch.  The  many  valued 
function  K  (jC)  Is  continuous  for  arg  ^  =0,  but  the  single¬ 

valued  principal  branch  Is  not.  Jones  used  the  principal 
branch  and  hence  his  values  show  this  discontinuity  along 
arg  C  =0.  Incidentally,  he  uses ’rectangular  coordinates 
which  hide  the  significance  of  the  restriction  on  arg  C  • 

Since  aerodynamlclsts  seem  to  be  most  interested  In 
the  region  around  arg  C  =0,  there  Is  something  to  be  said 
for  the  Luke  and  Dengler  approach.  That  Is,  their  extension 
of  the  definition  of  c(k)  is  the  only  one  which  agrees  with 
equation  B-1  for  Im(k)<  0  and  at  the  same  time  Is  the  analytic 
continuation  of  equation  B-1  across  the  half-line  arg  k  =  0. 

It  should  be  emphasized,  however,  that  this  approach  still 
leaves  one  with  a  discontinuity  along  the  negative  real  axis, 
l.e.  for  arg  C  =  ±  tt  (if  one  desires  a  single-valued  func¬ 
tion)  .  Since  Luke  and  Dengler  do  not  discuss  the  branches  of 
functions  Involved,  there  Is  some  ambiguity  In  their  defini¬ 
tions  . 

Moreover,  the  physical  meaning  of  equation  B-5  for 
Im  ( 0  >  0  Is  not  at  all  clear. 


Hence,  the  clrculatlonal  part  of  the  lift  force  caused  by  a 
unit  step  change  In  translational  velocity  Is  given  In 
general  by 


L(  s)  =  2p7TUb  C  j2f(  s)  +  J  C 


eKis  -  x)  (b-16) 


This  result  agrees  with  Blspllnghoff ,  etc.  (reference  B-7) 
for  a  lift  force. 
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